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I. INTRODUCTION 

The sign problem remains one of the biggest challenges 
of lattice field theory till this day. It is caused by a highly 
oscillatory and complex path integral measure after intro- 
ducing a chemical potential fi. This renders many the- 
ories, like four-dimensional Quantum Chromodynamics 
(QCD), inaccessible by common simulation algorithms in 
wide regions of the phase diagram. Many solutions have 
been proposed, see e.g. [1—7]. However, reliable numer- 
ical calculations in theories with a severe sign problem 
remain extremely challenging. 

One of the approaches to the sign problem is stochas- 
tic quantization [8], i.e. a complex Langevin evolution 
[9]. For a more detailed introduction see [10], and for an 
overview see [11]. It has been applied to many theories 
suffering from a severe sign problem, often successfully 
[12-18]. Furthermore it was applied outside the context 
of finite density calculations [19-21]. However, there are 
cases known where the complex Langevin evolution con- 
verges towards unphysical fixed points [22, 23]. Start- 
ing from early studies of complex Langevin evolutions 
[24-26], till this day the convergence properties of the 
generalization of Langevin equations to the complex case 
are not well-understood. In this work we focus on the 
stochastic quantization of the Thirring model at finite 
density in 2 + 1 dimensions, see also [27]. This serves as 
a toy model for the matter sector of QCD. Furthermore, 
the 2 + 1-dimensional model appears in effective theories 
of high temperature superconductors and graphene, see 
e.g. references given in [28]. 

Here we extend the + 1-dimensional studies carried 
out in [10]. In 2+1 dimensions we lose the analytic bench- 
marks which facilitated the interpretation of the numer- 
ical results with stochastic quantization in + 1 dimen- 



sions. Still, we can exploit the similarities to the lower- 
dimensional case. Further benchmarks of our numerical 
results are given by those in the heavy dense limit, as 
introduced in [29]. This limit describes the regime of 
large fermion masses and large chemical potentials. In 
addition we evaluate indirect indicators, namely the con- 
sistency requirements derived in [30] and the analyticity 
of the fermion condensate at /i 2 = 0. 

The paper is organized as follows: In section II we 
introduce a generalized Thirring model in the contin- 
uum and on the lattice. We also formulate the associ- 
ated Langevin equation. In section III we employ the 
heavy dense limit to derive numerous analytical results. 
At the end of the section we introduce indirect indicators 
of correct convergence, namely analyticity of observables 
at /i 2 = and a set of consistency conditions. In section 
IV we discuss numerical results and use the analytical re- 
sults among other indicators as a benchmark to evaluate 
the complex Langevin evolution. Finally in section V we 
discuss and summarize our findings. 

II. THE THIRRING MODEL 

A. Continuum formulation 

We begin with a short recapitulation of the model, 
which we introduced in [10]. It is a generalization of the 
historical 1 + 1 dimensional Thirring model [31], but for- 
mulated in d dimension with Nf fermion flavors at finite 
density. The Euclidean Lagrangian reads 

N f 
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Here i = 1, . . . , Nf is a flavor index, nn and /Lt« are the 
bare mass and bare chemical potential of the respective 
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flavor and g 2 is the bare coupling strength. In d = 2 + 
1 dimensions ^ and \& denote four-component spinors. 
While it is possible to use an irreducible 2-dimensional 
representation of the Dirac algebra, this representation 
does not allow for a chiral symmetry in the massless case. 
See [28] for details and a discussion of symmetries. The 
7-matrices satisfy {7^,7,,} = 28^1. 

We reformulate the generalized Thirring model using 
an auxiliary field and integrate out the fermionic degrees 
of freedom. We find for the partition function 



Z= j 3>A ^Qdet-fQJ e~ 

l/T 

S A = N f (3 J dt j d^xAl 



Sa 



<2>Ae 



(2) 



with temperature T, fermionic term Ki = $ + iA' + mi + 
|Ui7o and S e g = Sa — X^Trlogi-Q. The fermion deter- 
minant obeys 

d*AKi(ji) = [det Ki{-ii*)]*, (3) 

yielding in general a complex action. The sign prob- 
lem can be avoided by taking the absolute value of the 
fermion determinant. This corresponds to an isospin 
chemical potential and is also referred to as the phase- 
quenched case. We also note that observables become in- 
dependent of the chemical potential /1 up to some thresh- 
old n c in the zero temperature limit. This goes under the 
name of the Silver Blaze problem [32] . 



B. Lattice formulation 

We consider the generalized Thirring model in 2 + 1 
dimensions on a space-time lattice with N t time slices 
and N s lattice points in spatial direction. Furthermore 
we require that Nt is even [10]. The spatial volume and 
the space-time volume are denoted by 



V = Nf 



ft = iV t iV s 2 



(4) 



We use staggered fermions, where the number of flavors 
is given by J\f. The introduction of a chemical potential 
follows the prescription by Hasenfratz and Karsch [33]. 
In the following all dimensionful quantities are measured 
in appropriate powers of a, so that we only deal with 
dimcnsionlcss parameters. 

In three dimensions TV staggered fermion flavors cor- 
respond to Nf — 27V continuum flavors [27]. In principle 
the partition function for a single continuum flavor can 
be obtained by taking the square root of the fermion 
determinant, but it is still under debate if this rooting 
prescription is consistent [34]. 

The lattice action we use reads 

S = J2 X t (x) Ki (x, y) x, (y) + ^E A l (*) ■ ( 5 ) 



Here x% an d Xi denote staggered fermion fields. The sum 
extends over x, y = 1, . . . , ft and v = 0, 1, 2. The fermion 
matrix reads 



d-l 



Ki (x, y) = \Y. £ » [(! + iA - eM "°<W,i/ 

!/=0 

- (1 - xA v (y)) e-^ a 8 x .o,y] + m t S xy (6) 
with staggered phase factor 

s u ( X ) = (-i)^:=>* (?) 

and v denoting a unit vector in ^-direction, compare with 
[27, 35]. We impose periodic boundary conditions in spa- 
tial and antiperiodic boundary conditions in temporal 
direction. The lattice partition function is like in the 
continuum in (2) given by 



j \[dA v {x) ! IJdetiQ j e~ s \ 



(8) 



where Sa = ^2 X „ A 2 (x) . The central observables 
in our analysis are the fermion density, the fermion con- 
densate, the energy density and the phase factor of the 
fermion determinant. In the following, sums over flavor 
indices are not implied. The fermion density of flavor i 
is given by 



1 /dlogZ 



ft \ / \ 

I he fermion condensate follows from 
1 (d\ogZ\ 1 



(9) 



{XiXi 



ft V drtii 



and the energy density reads 
'd\ogZ s 



ft 



(TrK^) (10) 



dN t 



+ Hi (th) ■ 



(11) 



V,IH 



Usually we normalize the latter one to (si) (/x = 0) = 0. 

The phase factor of the fermion determinant is defined 
by det K = |det K\ exp (i</>) with phase <j>. For AT degen- 
erated staggered fermion flavors we set Ki = K. The 
expectation value of exp (ihf (j>) [36, 37] follows in the AT 
flavor phase-quenched theory from 



<^% q =f^[u,i] 



(12) 



where is given by (8) and Zff denotes the phase- 
quenched partition function. The expectation value of 
the fermion phase factor is a measure of the sign problem, 
where smaller values indicate a more severe problem. 
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C. Complex Langevin equation 



and 



To deal with the associated sign problem we apply 
stochastic quantization [8], namely a complex Langevin 
evolution [9], to the Thirring model at finite density. To 
this end we determine the stationary solution of the cor- 
responding Langevin equation, which reads 



(Re <q v 0, 6) Im rj CT (x', 6')) = 0. 



(20) 



d_ 



A 



v {x,Q)=- * S * [A l, +y/2T h ,(x,Q) 



8A V (x,Q) 



(13) 



Here denotes fictitious time and r\„ (ar, 0) is a Gaussian 
noise with 

fa (x,Q))=0, 
( Vl/ {x, O) rj a (x', 6')) = S va & (x - x 1 ) 8 (6 - 6') . (14) 

We use an adaptive first order stepsize integration scheme 
[27, 38, 39], where the stepsize e_t is adjusted with respect 
to the modulus of the drift term. 

For the numerical treatment of the Langevin equation 
we use a first order integration scheme. Our discretiza- 
tion of (13) reads 

A v (x,e + e L ) = 

A, (x, 6) + e L D v (x, 6) + V2^T] V (x, 6) , (15) 

where D v (x,Q) = —dScs/dA v (x,Q) is the drift term 
and ei, the (adaptive) integration stepsize. The drift term 
takes the explicit form 

D u (x,Q) = -M$A V (aj,0) 

+ ~e v {x)Y,Kr\x + v,x)e^ 



+ -e y (x + v) K% 1 (x,x + u) e 



— Mi<5„0 



(16) 



see also [27]. The stepsize is updated after each integra- 
tion step according to 



(L = £L (6) 



max Xjl , \D V (x, 6)| 



(17) 



Typically we use 8 = 1CP 2 , see also [10]. 

For a correctly converging complex Langevin evolution 
we can generalize the real noise to an imaginary noise 
in (13) while keeping expectation values unchanged, see 
also [27]. We parametrize the complex noise using the 
replacement 

r] v (x, 6) yJX + 1 Re r] u (x, 6) + i\/ZIm r\ v (x, 6) (18) 
with 1^0. Furthermore we require that 

(Re 77, (x, 6) Re 77, (x',6')) 

= (lmr] u (x, 6)Im?7 CT (x',@')) 

= S va S (x - x') 8 (& - 6') (19) 



While maintaining numerical stability, we check if ob- 
servables turn out to be independent of I. 



III. ANALYTICAL RESULTS 

In the following we derive approximate expressions 
for some observables in the lattice Thirring model in 
d dimensions. We begin with a hopping parameter ex- 
pansion and then take the so-called heavy dense limit 
in section IIIB. This renders the model effectively one- 
dimensional and allows us to obtain a simple expression 
for the fermion determinant. Later we discuss an ex- 
tended version in section IIIF. 



A. Hopping parameter expansion 

Applying a hopping parameter expansion to the 
fermion determinant in (8) yields 



det K 

rn ^ 



niK 1 - 1 ^) 

e. {c e } 



(21) 



Here n = 1/ (2m) is the hopping parameter, {Cg} are 
closed contours of perimeter £, n the number of times a 
given contour is traced out and Pc e is the product of 
hopping terms M (x, y) at /x = along the given contour 
Cg. The matrix elements M (x,y) read 



M (x, y) 



E 

u=0 



e v (x) [(1 + \A V (a;)) 8 x+c , iV 

-(l-iA^y))^.^], (22) 



compare with the fermion matrix in (6). Furthermore 
we represented the dependence on the chemical potential 
explicitly via 



7cv 



where 



•exp(//AT t )] 



6C e e Z 



6c e 



(23) 



(24) 



denotes the temporal winding number of the path Ce 
counted in positive direction. The minus sign in (23) 
stems from antiperiodic temporal boundary conditions. 



B. Heavy dense limit 

The heavy dense limit projects out the leading contri- 
butions of the hopping parameter expansion in the limit 
of a large mass m and large chemical potential /i. This 
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limit was introduced in [29] and employed e.g. in [15, 40- 
45]. It is denned by 



. — > 0, fi — y oo, ne^ fixed. 



(25) 



In this regime the fermion determinant is dominated by 
contributions from Polyakov loops in positive time di- 
rection with t)G = 1, see (24). Due to the absence of 
spatial paths, the model is effectively one dimensional. 
The fcrmionic contribution det /C in this limit reads 



det/C 



777 



n( i+ ^) 



(26) 



CEP 



where P denotes the set of Polyakov loops and £ = ( Nt 
with ( = Ke^. Let ^ x denote a Polyakov loop 



N t 



JI(l + L4 (t,x)) 



(27) 



starting and ending in the space point x. We can then 
express the fermion determinant in the regime of m and 
fi being large by 



det/C 



n 



i+^n( i+ii » (*> x )) 



(28) 



where O is the space-time volume defined in (4). Note 
that in this limit relation (3) is violated. Due to this 
approximation we will find (n) ^ =0 — > only for m — > oo. 
In section IIIF we introduce a modified version of this 
limit, which preserves this relation. 



C. The case of one flavor 

The simplified fermion determinant in (28) allows the 
derivation of an analytical expression for the partition 
function. We begin with TV = 1 flavors and later gener- 
alize to more flavors. In the heavy dense limit we find for 
the partition function 



(1 + 0' 
(2nf 



Hence the fermion density and condensate read 

' , {xx) 



1 



(29) 



(30) 



We point out that the exact + 1-dimensional results 
derived in [10] reproduce the above results in the corre- 
sponding limit. The normalized energy density reads 



{e = 



(C-K^logK- 1 
(1 + K 1/T) (! + £)• 



(31) 
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(a) Plot of the fermion density (n) 
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(b) Plot of the fermion condensate (xx) 

Figure 1: The phase structure in the heavy dense limit 



where we identified the temperature with T — N t . Fur- 
thermore we find 



£V C 

T 2 (1 + 0' 



> 0. 



(32) 



for the heat capacity. The corresponding mechanical 
equation of state reads 



P = log 



2k \p) 



(33) 



with P denoting pressure. 
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Figure 1 shows the phase structure in the heavy dense 
limit. Note that /3 dropped out in most considered ob- 
servables. For large N t we find two well-separated phases, 
where the system is in a condensed phase for large fj,. 

In the limit of vanishing temperature we find 



( n ) T =o = (^- Mc) , 
(xx) T =o = 2k0 (Mc - M) , 
(z)t=o = lo S (M - Me) ■ 



(34) 



where is the Heaviside step function and the critical 
chemical potential onset is found to be /x c = log (2m) . 
Note that [i c > in the heavy dense regime. The model 
clearly exhibits Silver Blaze behavior as mentioned in 
section I. 



D. The case of two flavors 

We continue with determining the partition function in 
(8) for TV = 2 flavors in the heavy dense limit. For non- 
degenerated flavors we label the parameters k, /i and £ 
with a flavor index i £ {1,2} and find 



Z, 



(1 + £i + £2 + £162 A) 



(4ki« 2 ) 



where we introduced 



1 



1 



(35) 



(36) 



for brevity. For simplicity we quote the observables only 
in the case of degenerated flavors and a common chemical 
potential, i.e. Ki = k, fa = fj, and £j = £. The generaliza- 
tion to non-degenerated flavors is straightforward. The 
total fermion density is given by 



2g(l + £A) 
1 + 2fJ + £ 2 A 



(37) 



with (n) — > 2 for n — > 00 (if /3 + 1 1/2), while the fermion 
condensate reads 



(xx) 



4k (1 + 



l + 2£ + £ 2 A' 

The mechanical equation of state takes the form 



P = log 



(l + 2g + $ 2 A 

4k 2 



l/T 



d/2' 



(38) 



(39) 



Like in + 1 dimensions [10] we find a plateau in the 
density, condensate and the energy density for TV = 2. 
They become visible for couplings of the order j3 as 1/2, 
see figure 7. For the special base of /3 — 1/2 the density 
never goes into full saturation for fj, — > 00, but is stuck 
on the plateau. These structures can be understood in 
the + 1-dimensional continuum case [10], see also [46]. 



Phase-quenched case. We can also consider the parti- 
tion function with a phase-quenched fermion determinant 
as mentioned in section II B. Then the partition function 
Z\ q in the case of H = 2 degenerated flavors reads 



1 



(2«) 



2Q 



(l + 2£ + £ 2 A pq ) 



N t d/2 



-1 V 



with 



A pq = 



1 



(40) 



(41) 



The observables are of the same form as in the full theory 
with A replaced by A pq . 

Phase-factor. The previous results allow us to derive 
an analytical expression for the expectation value of the 
phase- factor of the fermion determinant in the two flavor 
theory. It serves as a measure for the severity of the sign 
problem. We find 



2i0\pq 



<e 2i *) 
where we defined 

V = lo; 



;rpq 

^2 



-V/V 



l + 2g + C 2 A pq 
1 + 2^ + e 2 A 



(42) 



(43) 



As expected, the sign problem gets more severe for larger 
lattices. In the limit of large chemical potentials the ex- 
pectation value approaches 



lim / P 2i< n pq 



2/3-1 
2/3+1 



(44) 



We see that the choice of /3 has a significant impact on 
the severity of the sign problem. In the heavy dense limit 
we find that it is most severe for /3 = 1/2. 



E. The case of Af flavors 

Finally we generalize (29) to an arbitrary number of 
Af degenerated flavors. We obtain 



Zh = 



1 



E 

,fc=0 



k TT N t 



1-k 3 A/73 
2 ' 2' ~2~ 



(45) 



for the partition function in the heavy dense limit, where 
we introduced 



Afe = {ml 



(46) 
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d dimensions - Heavy dense limits - N t = 1 0, m = 0.5 
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Figure 2: The density in the strict and in the 
symmetric heavy dense limit 



Here T (z) denotes Euler's gamma function and U (a, b, z) 
is the confluent hypergcometric function of the second 
kind [47], also known as Rummer's function. We do 
not want to quote the observables here explicitly but 
as required they reproduce previous results for M = 1,2. 
From the expressions for the density and condensate we 
obtain the relation 

<*>_*(,_$). (47) 

In this general case we find up to J\f — 1 intermediate 
plateaus in the observables. 



F. Symmetric heavy dense limit 

The heavy dense limit given by (25) violates relation 
(3) . In the following we consider a model which restores 
the determinant symmetry, but is not a strict limit of the 
Thirring model, see [15]. Besides the Polyakov loops in 
positive time direction with their characteristic exp (/i)- 
dependence, we will also consider the ones in negative 
Polyakov direction with a exp (— /^-dependence. In this 
sense we speak of a symmetrized version of the heavy 
dense limit. 

In this model the fermion determinant det & has the 
form 

det£=n(l + it II (1 + i^o (x, t))^ 

xJJ^l + ^II(l-L4o(x,t))J. (48) 

where we introduced £/ = (/ce* 1 )^* and £j = («e~ M ) *. 
The respective partition function reads 

2f ym = ( % ) 2 (1 + it + it + A sym ) y (49) 



(2«r v p 



with shorthand notation 
A 



isym 



1 

1 + f 



N t 



(50) 



For small \i the observables now receive additional con- 
tributions. The density reads 



(») = 



io-it 



i+it + ib + A sym 



(51) 



and vanishes exactly at // = 0. The fermion condensate 
is given by 



(xx) 



2k (1- A sym ) 



i + it + ib + A sym 

and the equation of state takes the form 



P = loe 



(l + Q + 6 + A sym ) 
2k, 



1 ' /2vr^ 1 

J 



(52) 



(53) 



The quantitative differences between the ordinary and 
the symmetric heavy dense limit are small in the regime 
of large fermion masses and large chemical potentials due 
to the exponential suppression of the additional contribu- 
tions. For small m they alter the behavior of observables 
massively, see figure 2. However, for very small m the 
results again remain unphysical. Spatial contributions to 
the fermion determinant become important, thus render- 
ing the model invalid. 



G. Analyticity in /i 2 

For observables which are even in [i we can consider 
the analytic continuation to /i 2 < 0, corresponding to 
an imaginary chemical potential. An example for such 
an observable is the fermion condensate (Xx)- -Due to 
relation (3) the theory is free of a sign problem in this 
case and we can make use of a real Langevin evolution. 
For /i 2 > we employ a complex Langevin evolution. 

Assuming that the complex Langevin evolution is cor- 
rect, (xx) should be analytic in /i 2 . Then any non- 
analytical behavior at fi 2 = has to be caused by incor- 
rect convergence. We will check this criterion in section 
IV F, compare also with [13]. 



H. Consistency conditions 

In [30, 48, 49] the authors derived an infinite tower 
of identities able to indicate correctness of complex 
Langevin evolutions. They could show that for all entire 
holomorphic observables O the relation (L„ (x) O) = 
holds if and only if the evolution is converging correctly. 
Here 



L v (x) 



&A V (x) 



+ D v (x) 



&A V (x) 



(54) 
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is the Langevin operator. 1 A simple finite set of observ- 
ables is defined by O v (x, k) — exp (ikA u (x)) with k £ Z, 
yielding a necessary condition for correctness. In section 
IV F we will check if 

(L v O v ) = (ik [ik + D v (x)] e ikA ^) = (55) 

holds for a wide range of parameters. 

IV. NUMERICAL RESULTS 

A. Implementation 

We solve the Langevin equation numerically using the 
adaptive stepsize algorithm described in section II C. The 
implementation uses the GNU Scientific Library [50]. 
The code allows the calculation of the density, the con- 
densate, the phase-factor of the fermion determinant and 
the consistency conditions of section IV F. Furthermore 
there is the option to carry out calculations at purely 
imaginary chemical potential and in the phase-quenched 
case. 

To numerically determine observables we begin with a 
hot start, i.e. we randomly initialize the auxiliary field. 
This is followed by 5000 thermalization steps. To eval- 
uate the observables we sample the field configuration 
typically O (l0 4 ) times. Like in the + 1-dimensional 
case [10], we estimate the error with a bootstrap analysis 
[51] as the resulting error bounds generally proved to be 
more reliable. However, in some cases we observe that 
the actual error is underestimated. 
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B. Comparision to phase-quenched case 

In figure 3 we find a typical example for the density and 
the condensate in the full as well as the phase-quenched 
case. The error bars are included, but are often too small 
to be spotted with the naked eye. Although we only have 
N t = 8 time slices, we can already find hints for Silver 
Blaze behavior. The onset to the condensed phase for 
the full theory generally lies higher than in the phase- 
quenched case. 

The phase factor of the fermion determinant can serve 
as a measure for the severity of the sign problem and 
can be found in figure 4. We see that for large /3 it 
is less pronounced, but still quickly approaches zero for 
increasing fi. 

C. Evaluation in the heavy dense limit 

In the heavy dense limit the coupling factorizes, thus 
rendering the density and the condensate independent of 



(b) Condensate (xx) 

Figure 3: Full and phase-quenched numerical results for 
0=1,10 
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Figure 4: Severity of the sign problem 

In this section summation over v is not implied. 
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2 + 1 dimensions - Full theory - m = 100, N t = 8, N s = 8 
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Figure 5: Numerical results and heavy dense limit of 
the condensate for large masses 



(b) Fermion condensate (xx) 

Figure 6: Observables as function of (3 



D. Coupling parameter dependence 



(3. However, in figure 5 we can see that even for extremely 
large masses there is still a dependence on (3 in the numer- 
ical results obtained by a complex Langevin evolution. 
Furthermore, the gap between numerics and the heavy 
dense results seem to persist even when increasing m by 
several magnitudes. For smaller (3 this deviation is more 
pronounced, while for increasing (3 we can see how the 
curves are approaching each other asymptotically, similar 
to the + 1-dimcnsional case [10]. 



It is difficult to separate the different contributions to 
this gap as we are comparing an approximate result with 
an algorithm, which correctness we aim to check. In the 
+ 1 dimensional case we already found a real deviation 
to theoretical results for small /3, and something similar 
might also happen in the 2 + 1-dimensional case. 



In figure 6 we plotted the observables as a function of 
f3. For small /3 it is more difficult to obtain reliable re- 
sults due to spikes caused by numerical instabilities. For 
very large (3 the results from the full theory approach the 
phase-quenched calculations. As the heavy dense limit 
makes no predictions about the f3 dependence, it is diffi- 
cult to compare with analytical results. 



E. Several flavors 

in figure 7 we can find the density and the conden- 
sate in the case of M — 2 staggered fermion flavors at 
a coupling of f3 — 0.6. As we saw in section HID and 
HIE these observables show up to M — 1 plateaus in a 
certain range of j3. However, we observe that the numer- 
ical results obtained by a complex Langevin evolution do 
not reproduce these plateau structures at any value of j3. 
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2 + 1 dimensions - N f = 2, p = 0.6, m = 3, N t = 8, N s = 8 



2 + 1 dimensions - Full theory - p = 3, m = 1 , N t = 8, N s = 




Chemical potential n 
(a) Fermion density (n) 




Chemical potential fi 



Figure 8: Condensate (xx) with imaginary noise X 



2 + 1 dimensions - N ( = 2, p = 0.6, m = 3, N t = 8, N s = 8 




for X = 0, TV = ft = m = fi 
condition. On the other hand 



1 for a violated 



2 3 
Chemical potential u 



(b) Fermion condensate (xx) 

Figure 7: Observables for J\f = 2 degenerated flavors 



Only when considering very large /3, the theoretical curve 
is free of these structures and the complex Langevin evo- 
lution approaches the correct result. This finding is in 
agreement with [10]. 



F. Consistency conditions 

The check of the consistency conditions of section III H 
for the O v (x, k) observable yields several violated condi- 
tions. Without loss of generality we restrict ourselves 
to the case of x 1 = (1,1,1). For improved statistics we 
consider a 6 3 lattice, where we sampled the field con- 
figuration 5 • 10 times. As an example we quote here 
excmplarily 



Re (L Q (x) Oo (x, k)) 
Re (Li {x)Oi {x,k)) 



0.2258 ± 0.0192, (56) 
-0.0244 ± 0.0049 (57) 



Re (L (x) Co (x, k)) = 0.0080 ± 0.0058, (58) 
Re (Li (a;) 1 (x, k)) = -0.0004 ± 0.0059 (59) 

for X = 0, TV = m = ^ = k = 1, /? = 10 is an example for 
a condition which is compatible with a vanishing value. 
As usual denotes the number of staggered fermion 
flavors, X the imaginary noise, /3 the inverse coupling 
constant, m the mass and \x the chemical potential. 

If we take the error bounds estimated from the boot- 
strap analysis seriously, we observe several violations of 
the consistency conditions in 2 + 1 dimensions. This is in 
agreement with the + 1-dimcnsional case, where many 
consistency conditions seem to be not fulfilled. We take 
this as a hint that the complex Langevin evolution in this 
simple implementation might not converge correctly. 



G. Imaginary noise 

We discussed in section II C how we can generalize 
from a real to an imaginary noise term, assuming cor- 
rect convergence of the complex Langevin evolution. In 
this case observables should turn out to be independent 
of the noise term X. However, in figure 8 we can find clear 
evidence for a dependence. Furthermore we observe for 
almost all I > a significant larger deviation to analyti- 
cal results. Hence it was also not possible to fine-tune X 
to a value so that the complex Langevin evolution yields 
correct results. 



H. Analyticity in fi 

In figure 9 we find the condensate (xx) as a function 
of fi 2 . While for /i 2 > the numerical evaluation em- 
ploys a complex Langevin evolution, for /i 2 < we use 
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2 + 1 dimensions - Full theory - m = 1 , N t = 8, N s = 8 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 

Chemical potential fi 2 



(a) Mass m = 1 at varying /3 



2 + 1 dimensions - Full theory - p = 3, m = 0.2, N t = 8, N s = 8 
0.25 | . , , , . 1 . 1 , 1 . 1 , 1 — 




-4-3-2-101234 
Chemical potential u 2 



(b) Mass m = 0.2 at /3 = 3 

Figure 9: Condensate (xx) as function of fi 2 

a real one. For small ft the numerical evaluation tends 
to be unstable. Like in the + 1-dimensional case (xx) 



is analytic within the numerical accuracy. This suggests 
that the method can work sufficiently well for small fi. 



V. CONCLUSIONS 

In this letter we have extended our + 1-dimensional 
investigation of the generalized Thirring model in [10] to 
the 2+ 1-dimensional case. Our numerical results are sim- 
ilar to those obtained in + 1 dimensions. The complex 
Langevin evolution allows us to evaluate observables in 
a straightforward way, despite the severe sign problem. 
However, our investigation suggests that in some cases 
it does not converge towards the physical theory. This 
was indicated by several violated consistency conditions, 
a gap to analytical predictions in the heavy dense limit 
and the absence of plateaus in the observables for JV > 1 
flavors. 

Despite these observations, we found the condensate 
(Xx) to be analytical at y? = 0, suggesting that a com- 
plex Langevin evolution seems to work for small fi. Also 
for large ft the heavy dense limit results could be re- 
produced. In these regimes we then have an appealing 
method to tackle the sign problem in the generalized 
Thirring model. 

Further investigations have to deal with the question 
of how to address the aforementioned problems. In par- 
ticular, coordinate transformations as suggested in [52] 
and gauge cooling procedures like the one employed in 
[53] might allow a stabilization of the complex Langevin 
evolution. 
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